Landauer conductance of tunnel junctions: Strong impact from boundary conditions 
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We present model calculations for the Landauer conductance of tunnel junctions. The tunneling 
of free electrons through a rectangular potential barrier is considered. The conductance of a finite 
number of barriers was calculated using a transfer matrix method. The periodic arrangement of the 
same barriers was described by a Kronig-Penney model to calculate the band structure and from that 
the conductance of a point contact in the ballistic limit. The comparison of both results elucidated 
the importance of the boundary conditions. Caused by resonant scattering in the superlattice the 
conductance is overestimated by the order of 1/t, the transmission coefficient of the single barrier. 
In the case of metallic multilayers these interferences are of minor importance. In conclusion, the 
application of the Landauer formula to periodic lattices to describe the tunneling conductance of a 
single barrier is not correct. 

PACS numbers: 73.23.Ad,73.63.Hs 



I. INTRODUCTION 

The discovery of Giant Magnetoresistancei*^ in 1988 
initiated a renewed interest in the phenomenon of spin 
dependent tunneling. First Julliere in 197E 3 succeeded 
to measure the spin polarized current in a tunnel junc- 
tion composed of two ferromagnetic layers separated 
by an insulating barrier. Using improved experimen- 
tal techniques it is now possible to produce high qual- 
ity tunnel junctions showing a remarkable tunneling 
magnetoresistance^*Si2i. A widely used model to inter- 
pret the results in terms of the polarization of the ferro- 
magnetic leads was proposed by Jullierei It was based 
on earlier studies of Meservey and Tedrow on tunneling 
between superconductors and ferromagnets&. The model 
assumed that the tunnel current is proportional to the 
product of the effective spin-dependent tunneling density 
for both ferromagnetic leads. Using parabolic band mod- 
els the influence of the barrier height and width could be 
included in the theoretical investigation al' 1 ? . 

The electronic structure of the magnetic electrodes is, 
however, quite different from the free-electron picture. 
This could be incorporated in the considerations using 
tight-binding model hamiltoniansiii^. Provided by the 
simplicity and less demanding numerical effort of these 
models the influence of defects on the tunneling current 
could be considered, too. The up to now most sophisti- 
cated models use the self consistent potential and elec- 
tronic eigenstates of the whole tunnel junction to de- 
termine the transport propertiesiiiiii. A well estab- 
lished formalism to treat tunneling conductance follows 
the early ideas of Landauer and Biittike»i£ii£. The cru- 
cial property to determine is the tunneling probability 
of the reservoir states through the barrier. It is mostly 
calculated using the Kubo-Greenwood approach based 
on Green's functions or a multiple scattering formalism 
based on the electronic eigenstates. In the experiments 
the thickness of the electrodes is much larger than the 
electron mean free path representing homogenous reser- 



voirs and justifies the assumption made by Landauers 
formalism. Assuming semi-infinite boundary conditions 
for the electrode hamiltonian the theoretical transport 
description is very close to the experimental conditions. 
Density functional theory based self-consistent electronic 
structure codes assume sometimes a three-dimensional 
periodicity of the system. This means that the tunnel 
junction geometry is replaced by a periodic system con- 
sisting of electrodes of finite thickness separated by tun- 
nel barriers with thicknesses used in the experiment. The 
implications for the self consistent potential and charge 
densities are small due to the small screening length 
in the metallic regions. Considering metallic multilay- 
ers Landauers formula was applied to the Co/Cu(001) 
systemiSiiS. In this case the theory describes the current 
through a (peculiar) point contact made of a thin pillar 
and electrodes both with the same periodic superstruc- 
ture in current direction. 

This paper is aimed at the importance of the bound- 
ary conditions for the description of transport proper- 
ties of tunnel junctions. The initial idea was to deter- 
mine starting from the eigenstates of a periodic array of 
junctions the current through one of them in a geome- 
try with infinite leads (electrodes) on both sides. The 
forthright attempt would be to adapt Landauers formula 
to the periodic lattice geometry assuming perfect trans- 
mission for all Bloch states^. This attempt assumes 
a perfect contact of the barrier to reservoirs with elec- 
tronic states which are caused by the periodic lattice. 
The conductance is determined by the contact resistance 
only 27 . One could suppose that the superlattice quanti- 
zation caused by the reflection at the barriers should com- 
prise the transmission probability occurring in the Lan- 
dauer formula for a single barrier. By this assumption 
the determination of the transmission probability of the 
(semi-infinite) electrode eigenstates would be avoided. 
The discussion later on will show the limitations of this 
approach and the failure in the case of tunnel junctions. 
To determine the conductance in the tunneling regime 
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FIG. 1: Geometrical arrangement of rectangular barriers of 
height U and width L, a) single barrier case, b) finite number 
N of barriers (shown for N = 2) with a center-center distance 
of a, and c) periodic lattice 



barriers of width L with a constant distance a , and zero 
potential in the valley regions between the barriers and 
on both sides as shown in Fig^>). The barrier height 
U = .38Ry and width L — 16.9a.u. were fixed, if not oth- 
erwise stated. In the periodic lattice the arrangement of 
barriers is translational invariant in z-direction. Empha- 
sizing the periodicity of a system throughout this paper 
we refer to the configuration in current direction which 
is the z-axis. For the directions perpendicular to the cur- 
rent a constant potential is assumed. So, the momentum 
perpendicular to the z-axis k|| is a good quantum number 
and an effective one-dimensional problem is to be solved 
for every k|| . 

First, we consider a finite number N of barriers on both 
sides in contact with free electron reservoirs with zero 
potential and equal Fermi level E F . The conductance 
in zero bias limit is obtained using Landauer's formula^ 
adding up all states at the Fermi level 



G(N) = T £ rW(^-A|| 2 ), (1) 



the determination of the transmission probability is im- 
peratively necessary. To extract these quantities from 
a periodic lattice calculation a so-called inverse Kronig- 
Penney-model was analyzed. The reversibility of them 
could be to demonstrated obtaining the scattering prop- 
erties of one barrier from the eigenstates of the periodic 
system by means of a transfer matrix formalism^. 

Using a parabolic band model a piece wise constant po- 
tential in z-direction was assumed. The transition from 
a metallic to an insulating junction was studied by ad- 
justing the Fermi level with respect to the barrier height. 
Similar free electron models have been used by several au- 
thors to describe successfully the transport properties of 
magnetic multilayersSi* 2 ^ 2 ^* 2 ^ 2 ^. Starting from one rect- 
angular barrier we will evaluate the transmission proba- 
bility and the resulting conductance of the free-electron- 
like reservoir states through a finite number N of barri- 
ers. In the limit of large numbers N a coincidence with 
the properties of the periodic system could be expected, 
but could not be approved. We will discuss the implica- 
tions of the periodic boundary conditions on the trans- 
port properties and the limitations of the periodic lattice 
approach. The assumption about the shape of the poten- 
tial are not important at all for the general conclusions 
about the applicability of Landauer's formula. 



II. MODEL AND TRANSPORT EQUATIONS 

To emphasize the influence of the boundary conditions 
we have chosen a model system for the tunneling geome- 
try. The 3-dimensional potential was chosen to be piece- 
wise constant. The potential inside the barrier was set 
to U, outside it was set to zero. The modulation of the 
potential is along the z-direction. The thickness of the 
barrier was fixed to L. The single barrier case is sketched 
in Fig^,). A finite stack of barriers was simulated by N 
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%- is the conductance quantum, T^ N ^(E F — k^ 2 ) denotes 
the transmission probability of state k|| with an energy 
equal E P derived from the complex transmission coeffi- 
cient fl N \ Due to the simple potential shape the trans- 
mission coefficient depends on the kinetic energy in z- 
direction only. It was assumed that both reservoirs are 
identical and a single channel exists for a given fcj| . The 
factor 2 for spin degeneracy will be omitted throughout 
this letter. The neglection of the reflection probability 
with respect to the primary formulation of Landauer cor- 
responds to the conductance value obtained by a two 
point measuremen t 17l26 i. The difference to four point 
probe experiments arises from the contact resistance 27 
and is negligible if the transmission T is much smaller 
than unity. In the metallic regime for T near to unity 
this difference can outweigh the calculated part of the 
conductance. Our considerations are mainly focused on 
the tunneling regime. 

Setting A the cross section of the point contact the 
sum in eq. can be transformed 



r (N) = e A 
h4n 2 



f d 2 k||TW(B P -fc|| 2 ) (3) 

Jk« 2 <E F 



which can be written as an energy integral using the rota- 
tional symmetry of k|| <i 2 k|| 



q(n) = 



h 47T 



27rc?fc||fc|| and E 
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Ec 
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It is Evident, that the behavior of the conductance 
in dependence of the filling of the reservoirs, given by 
irrelative to the barrier height U , is simply related 
to the transmission probability T^ N \E) derived from 
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a one-dimensional model. The transmission coefficient 
& N \E) will be calculated using a transfer matrix for- 
malism sketched in the next paragraph. 

For a periodic arrangement of barriers we assume a 
perfect transmission for all Bloch states. That is, we 
consider now a periodic crystal described by a prolonged 
unit cell in z-direction containing one barrier and one 
well region (super cell approach with model potential). 
In the description of Landauer this corresponds to a bar- 
rier in contact with reservoirs which have a modulated 
potential as well. In this case the boundary conditions 
at both sides of the barrier are not free electron like, and 
the states in the barrier region have to match the eigen- 
states of the modulated reservoir potential. Many ab 
initio electronic structure codes are adapted to a lattice 
symmetry with 3-dimensional periodicity. So it seems 
to be desirable to use these eigenstates to calculate the 
transport properties even in the case of tunneling bar- 
riers. The conductance is given by an integral over the 
3-dimensional Fermi surfac o 18 i 19 
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It is a sum enclosing the projections S v of the Fermi sur- 
face sheets v into the plane perpendicular to the trans- 
port direction n (z-axis) 
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To compute this integral for our model system the band 
structure has to be evaluated at the Fermi level. A 
Kronig-Penney-model from the textbooks was used and 
is sketched briefly in the next paragraph. 



III. SCATTERING PROPERTIES AND FERMI 
SURFACE 

The derivation of the transmission coefficients of a fi- 
nite number of barriers should be sketched briefly. We 
start with a single barrier, and assume one incoming and 
one outgoing solution on both sides of a symmetric bar- 
rier (Fig. \^). The wave function amplitudes on both 
sides entering the general solution of the Schrodinger 
equation are connected by a 2x2 transfer matri:s2& 



M = 



1/t f*/t* 
f/t 1/t* 



(7) 



with the complex transmission and reflection coefficients 
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with E = k 2 and U — E = n 2 and t(E) the absolute value 
of the transmission coefficient. The parameterization of 



M according to eq. J7J) is based on the assumption of a 
symmetric barrier {V(z) = V(—z)). 

The total transmission coefficient of a stack of N bar- 
riers can be calculated by successive application 
of the corresponding transfer matrices. The transfer ma- 
trix M of the single barrier was derived for a symmetric 
barrier at z = 0. The transfer matrices of the other 
barriers in the stack have to be transformed according 
to their shift in z-direction by appropriate D matrices. 
The transmission coefficient t^ N ' depending on E — k 2 is 
obtained from 

M (N) = D (N-l)/2^ MD) (N-l) MD (N-l)/2 

and 
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The behavior of i<W (E) is shown in Fig.0 t 1 (E) = t(E) 
describes the single barrier case. For N larger than one 
resonances occur for energies below U corresponding to 
bound states in the valley regions between the barri- 
ers. At the resonance the transmission coefficient reaches 
unity in line with the assumption of elastic scattering 
only. The maxima of the peaks are not covered by the 
chosen energy grid due to their tiny width. The impli- 
cations of additional inelastic scattering are discussed in 
the Iiterature22iiffli2i. In the lower panel of Fig.|2the be- 
havior of the transmission next to a resonance is shown 
on a finer energy mesh. The number of sub-resonance 
peaks within one resonance interval increases with the 
number of considered barriers. To obtain G^ N \E F ) an 
integration over t^ N '(E) from to E F has to be per- 
formed according to eq. @ . The main contributions are 
given by the resonance peaks which are a tiny fraction 
of the total interval (typical 10~ 5 to 10~ 3 ). An adaptive 
mesh refinement algorithm similar to Ref^ was devel- 
oped to improve the Simpson integration scheme up to 
the desired accuracy. 

To obtain the Bloch states of a periodic arrangement of 
barriers, sketched in Fig.^O, the Kronig-Penney model 
was used. The energy bands are obtained from an im- 
plicit equation for the perpendicular wave vector k of the 
superlattice 



cos(fca + 8) 
t 



(11) 



which depends on the perpendicular wave vector k|| by 
E = k 2 = E F — fc|| 2 , the transmission coefficient t = 
t(E), and phase shift 6 = 5(E) of a single barrier defined 
in eq. ©■ The set of fen values providing a solution of 
eq. (|ll|l defines the projection of the Fermi surface sheets 
required in eq. © 
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FIG. 2: Transmission coefficient v- N '(E) for stacks of TV rect- 
angular barriers (U = .38Ry, L — 16.9a.u.,a = 54a.u.), lower 
panel: behavior of t^(E) at resonance 



k 2 v< and k?, v> are the inner and outer radius of the annu- 

lus and are given by the condition cos ( k ^+ s ) — ±1. First 
the zeros and local maxima of cos(fca + <5) are determined 
for k S (0, \/E F ), and the interval was partitioned in a 
set of smaller ones. Every smaller interval providing one 
solution for eq. was treated by a scheme of nested 
intervals. 



IV. RESULTS AND DISCUSSION 

The behavior of t^ N \E) is shown for several num- 
bers N of identical barriers in Fig. [2] The upper panel 
shows the over-all dependence including the basic solu- 
tion t(E) = tW(E) for one barrier. For energies larger 
than the barrier height U a nearly perfect transmission 
is obtained with narrow anti-resonances with a reduced 
transmission. For energies lower than the barrier height 
U the general trend of t (N) (E) is given by [t(E)] N . Reso- 
nances with perfect transmission occur at energies related 
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FIG. 3: Landauer conductance G'"^ (E f ): a) general trend, 
the curves for N = 2 and N = 4 nearly coincide, b) G (2) {E F ) 
for different barrier distances a with higher energy resolution. 



to bound states in the valleys between the barriers, see 
the lower panel of Fig. [21 For N — 2 this corresponds to 
the results for double tunneling diodes without inelastic 
scattering 30 - 31 - 33 - 34 . Expanding t^{E) in the vicinity of 
the resonance energy E r a Lorentzian peak is obtained 



t {2 \E) = t 2 (E) 



2i(ka+S{E)) 



r 2 (E) 



(e - E r y + r 2 



(13) 



assuming a slowly varying 8(E) and t — t(E) <C 1 
in the vicinity of E r . E r results from the condition 
cos(k(E)a + 6(E)) = 0. The obtained value for T = 
t 2 k/(ra) is comparable to the resonance width induced 
by defect states in the middle of the barrier in the case 
a — L <C a 34 . For N > 2 barriers N— 1 peaks with perfect 
transmission appear around E r . The width of the total 
resonance is about 



A r = Ay/E~t/a 



(14) 



which is of the order < _1 larger than T and corresponds 
to the band width of the Bloch states in the periodic 
system. 

Fig. [31 shows the dependence of G^ N ' as a function of 
the Fermi level E F . The numerical results are given in 
units of the prefactor %- in eq. [3 The height of the 
barriers was fixed and small E F describes transport in 
the tunneling regime, but large E F with respect to U 
models the (free electron like) metallic regime with pe- 
riodic potential perturbations. For small band fillings 
(E F <§; U) the conductance increases step- wise caused 
by the sharp resonance structure of the transmission co- 
efficient tW{E) for N > 1. With increasing Fermi level 
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the resonance in appears in k-space first at the 

fc|| = point because the effective energy of the one- 
dimensional tunneling problem E = E F — k\ 2 is the 
largest for these states. The spacing of the transmission 
resonances with respect to energy scales with the inverse 
of the distance of the barriers a (see eq. (tT^t 'l. As can 
be seen in Fig. I^Jd) for larger barrier distances the tun- 
neling conductance of TV > 1 follows exactly the trend of 
G^(E F ) apart fr om a factor of the order 2. 

The conductance contribution being added by one res- 
onance in the integral in eq. Q can be estimated in the 
case of small t = t^(E r ) by the product of the aver- 
aged transmission probability tW and the width of the 
resonance interval. From the analytical expression for 

t( N )(E) (for t {2 \E) see eq. l|T3)l) the mean value of 
averaged over one resonance interval is nearly equal \t 
independent on N . The total width of the resonance in- 
terval is given by Eq. 1141 Consequently the contribution 
of one resonance is given by 



AG 
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— 1 = 

a 4 



E 



(15) 



This is exactly half of the value obtained by the conduc- 
tance integral for one barrier in eq. J3J estimated by the 
resonance spacing 



AE = 2Akk = 2-\[~E 



(16) 



times the transmission probability t 2 (E) . The resonance 
spacing is derived from the condition that the resonances 
appear at energies of bound states in the valley region be- 
tween the barriers assuming hard walls (Aha = tt). This 
estimation is independent on the number N of barriers 
and causes the same qualitative behavior of the conduc- 
tance for a chosen barrier distance a, nevertheless the dis- 
tance of the barriers a determines the height and width 
of the resonance steps in G^ N \E F ), compare Fig. Et>)- 
The cases N = 2 and N = 4 coincide exactly in Fig.|3Ji. 

In the following attention is drawn to eigenstates, 
Fermi surfaces and derived conductance integrals for the 
periodic lattice containing the above described model 
tunneling barriers. The Fermi surfaces of the periodic 
lattices are shown in Fig. For reasons of comparabil- 
ity the Fermi level was fixed and the barrier height was 
chosen to be 0.38, 0.76, and 1.52 times E F represent- 
ing the cases of a large, nearly equal and small Fermi 
energy E F with respect to U , respectively. The upper 
row in Fig. ^ shows the dependence of the single barrier 
transmission T« = t 2 {E F - fcjj) on fc||. The lower row 
shows the respective Fermi surfaces in a projection on 
the plane perpendicular to the current direction. The 
radius of the outermost Fermi surface sheet corresponds 
to the chosen Fermi level. The area of the Fermi surface 
sheets projection is a direct measure for the conductance 
contribution in eq. 10 • The strong relation of transmis- 
sion probability and the width of the Fermi surface sheets 
is evident. For nearly perfect transmission free electron 




FIG. 4: Transmission probability T(fc||) in dependence on ku 
(upper row) and supercell Fermi surface (lower row), E F = 
.25Ry, barrier height U is set .38 (a), .76 (b), and 1.52 (c) 
times E F 
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FIG. 5: Probability amplitude profiles along the z-axis of 



typical eigenstates of the periodic lattice, a) E F — fey 
b) E F - fe,| 2 « U, and c) E F ~ fe,, 2 < U 



> u, 



like sheets with tiny gaps are obtained (left panel, cen- 
tral part). For small transmissions narrow sheets occur 
(right panel). The width in direction of /c|| can be esti- 
mated from eq. (|11|) by linearizing the cosine, assuming 
a slow variation of t(E) and 8(E), and the condition that 
the absolute value of the r.h.s. is smaller than one 



A, 2t 

Ak= — 

a 



(17) 



It corresponds to the energy width A r of the transmission 
resonance mentioned in eq. (f 141) above. The center of 
the band is given by the condition cos(fca + 5) = with 
k 2 = E = E F — k ]{ 2 . 

The scattering properties of the single barrier are re- 
flected by the eigenstates of the periodic system and can 
be determined using an 'inverse' Kronig-Penney-modelSi. 
Typical eigenstates at the Fermi level for different kn are 
characterized in Fig. [S] by their probability amplitude 
profiles along the z-axis. E F was chosen to be larger then 
U (Fig.0] left panel). Fig. [S^,) shows an eigenstate close 
to the center of the Brillouin zone (fc|, «0, E F -k l{ 2 > U). 
The probability amplitude is slightly modulated by the 
periodic potential, because the transmission of the bar- 
rier is nearly perfect (t w 1). For states with a larger kn 
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FIG. 6: Conductance G in dependence on band filling 
of the reservoirs (E F ), for comparison the exact result for a 
single barrier G' 1 ' is given, (U = .38Ry, L = 16.9a.u.) 



the energy E F — fc|| 2 decreases and Fig. Et 1 ) shows a state 
close to E F — fc|| 2 = U. A resonance inside the barrier is 
a typical feature for these states with a moderate trans- 
mission coefficient (t < 1). States with a small effective 
energy E = E F — fc|| 2 <C U which corresponds to large 
k\\ are confined to the valley region between the barriers, 
caused by the strong reflection (t -C 1, see Fig. Eh). 

Figure shows the main result of these investigations, 
the dependence of the conductance G sc on the Fermi 
level applying eq. ©■ For comparison the exact result 
for the one barrier conductance G^ 1 -* and the transmission 
coefficient t(E) are given. G sc is in good agreement with 
G^ in the metallic regime (E F /U > 1). The restrictions 
caused by the neglection of the reflection coefficients R in 
the Landauer formula discussed above would rescale both 
curves in the same manner 26 . In the tunneling regime 
(E F <C U) the supercell conductance G sc overestimates 
the correct value by the order of inverse t, as it is visible 
in the logarithmic plot. 

The discrepancies can be understood comparing the 
corresponding contributions to G^ in eq. J3J and to 
G sc in eq. ijfjj for one resonance interval and the cor- 
responding Fermi surface sheet, respectively. One has 
to assume large a, so that t(E) and 5(E) do not change 
significantly in a resonance interval. A similar estima- 
tion as given above for AG^ in eq. H15|) can be done 
for G sc evaluating the contribution for a single Fermi 
surface sheet. Disregarding the prefactors of the integral 
in eq. 10 the contribution of one Fermi surface sheet is 
given by the circumference 2nk — 2n\f~E times the width 
Ak='f from eq. l(T7jl 

S v = 2ttVe — . (18) 
a 

The corresponding contribution AG^ in eq. 10} is given 



using eq. (jTTj|) by 

AG (1) AG (1) = AEt 2 = 2n/a^/Et 2 . (19) 

Comparing these two expressions including the prefac- 
tors in eq. J3J and eq. JHJ) the ratio of G sc and G^ can 
be estimated by - 2 r. The main contribution to the intc- 
grals is given by the Fermi surface sheet and transmission 
resonance interval, respectively, with the largest effective 
energy E F — k\ 2 for small fcii. So the ratio G sc /GW can 
be estimated by the transmission coefficient t(E F ) at the 
Fermi level. A closer look to the curves in Fig. El reveals 
the prefactor 2/irt. 

Another picture to illustrate the differences of both ap- 
proaches is the following. Using the supercell approach a 
perfect transmission for all periodic lattice Bloch states is 
assumed. In terms of the transmission coefficient t(E) in 
dependence on energy this means to set v- N > equal unity 
inside the resonance intervals for large N. (see Fig. Eb), 
and to zero outside the resonance regions. The resonance 
region for large N give exactly the position of the periodic 
lattice bands. The differences outside the resonances are 
negligible but the contribution of the resonances will be 
overestimated by far, as given by t( N h One has to keep 
in mind the logarithmic scale in Fig. Et>) and the small 
width of the N — 1 sub-resonance peaks. 



V. CONCLUSIONS 

Based on calculations of the Landauer conductance for 
tunneling barriers in different geometrical configurations 
the influence of the boundary conditions could be elu- 
cidated. Comparing the conductance of one and a fi- 
nite number of identical rectangular barriers we found 
the same qualitative behavior in the dependence on the 
band filling of the electrodes (Fermi level) with respect to 
the barrier height U. For Fermi levels smaller than the 
barrier height (tunneling regime) both results coincide 
despite a factor of 2 even for large numbers N. Using a 
supercell approach and calculating the conductance from 
the periodic lattice eigenstates assuming perfect trans- 
mission, one overestimates the exact result in the tunnel- 
ing regime by order t , the transmission coefficient of 
the single barrier for states at kn = 0. In the tunneling 
regime the small transmission and the resulting interfer- 
ences of the states in the reservoirs prevent a description 
of the single barrier conductance by a periodic lattice. 
The inclusion of the correct semi-infinite boundary con- 
ditions for the eigenstates of the reservoirs is crucial. 

In the metallic regime (E F > U) both results coin- 
cide and the approach was already successfully applied 
to metallic multilayer point contacts in CPP geometry^. 
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ics. P.Z. likes to acknowledge financial support by the 
Swedish Foundation for Strategic Research and fruitful 
discussions with Susanne Mirbt. 
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